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Using the continuous-time random walk (CTRW) approach, we study 
the phenomenon of relaxation of two-state systems whose elements evolve 
according to a dichotomous process. Two characteristics of relaxation, the 
probability density function of the waiting times difference and the re¬ 
laxation law, are of our particular interest. For systems characterized by 
Erlang distributions of waiting times, we consider different regimes of relax¬ 
ation and show that, under certain conditions, the relaxation process can 
be non-monotonic. By studying the asymptotic behavior of the relaxation 
process, we demonstrate that heavy and superheavy tails of waiting time 
distributions correspond to slow and superslow relaxation, respectively. 

PACS numbers: 05.40.Fb, 02.50.Ey, 76.20.+q 


1. Introduction 

The relaxation processes describing the transition of macroscopic sys¬ 
tems from one equilibrium state to another are the subject of great interest. 
This is mainly because the characteristics of these processes contain im¬ 
portant information about the mechanisms of relaxation. These relaxation 
processes are usually studied under the condition that a constant general¬ 
ized force, which has been applied for a long time, is abruptly switched off 
(for details see Ref. [T] ). In this case, one of the main characteristics of 
the relaxation process is an extensive thermodynamic variable conjugate to 
the generalized force. In particular, for dielectric relaxation (see, for ex¬ 
ample, ©Oil] and references therein) the pair, the generalized force and 
its conjugate variable, consists of the external electric field and the electric 
dipole moment of the relaxing system. Similarly, for magnetic relaxation 
m 0 in the corresponding pair consists of the external magnetic field and 
the magnetic moment of the system. 
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An important characteristic of any relaxation process is the relaxation 
law, i.e., the properly normalized (dimensionless) conjugate variable as a 
function of time. In the simplest case, when the rate of change of this 
function is proportional to its magnitude, the relaxation law is exponential. 
However, many systems exhibit anomalous, non-exponential, relaxation [8]. 
For instance, the slow magnetic relaxation was discovered in systems of 
single-molecule mmm and single-chain mmm magnets. 

There is a wide class of systems whose relaxation properties, including 
anomalous ones, are completely characterized by the individual properties 
of their structural elements (such as single-molecule magnets, single-chain 
magnets, single-domain ferromagnetic particles, etc.). In particular, this 
happens when the state parameter of each element evolves according to a 
dichotomous random process (see below for a more detailed discussion). In 
this paper, we derive a number of analytical results describing the relaxation 
of these so-called two-state systems. 

The paper is organized in the following way. In Sec. El we introduce the 
dichotomous model of relaxation of two-state systems, establish its connec¬ 
tion with the CTRW model, find the Fourier-Laplace representation of the 
probability density function of the waiting times difference and the Laplace 
representation of the relaxation law. The laws of biased and unbiased re¬ 
laxation of these systems characterized by Erlang distributions of waiting 
times are derived in Sec. [3j Here we show that under certain conditions the 
relaxation process can be non-monotonic. In the same section, by studying 
the asymptotic behavior of relaxation processes in two-state systems, we 
demonstrate that heavy and superheavy tails of waiting time distributions 
are responsible for slow and superslow relaxation, respectively. Finally, our 
findings are summarized in Sec. |H 

2. Dichotomous model of relaxation 

2.1. Model description 

We consider a relaxing system as a set of identical objects, each of which 
can be in one of two states randomly changing in time. This model is widely 
used to describe the relaxation processes in physical systems whose struc¬ 
tural elements are approximately characterized by two equilibrium states. 
Magnetic systems consisting of uniaxial single-domain ferromagnetic par¬ 
ticles represent an important class of such systems. The magnetization of 
each particle has two equilibrium directions, but due to the thermal fluctu¬ 
ations, its instantaneous direction can be arbitrary. In this case, the mag¬ 
netic relaxation is well described by using the stochastic Landau-Lifshitz 
or Landau-Lifshitz-Gilbert equation for magnetization dynamics and the 
corresponding Fokker-Planck equations for the probability density of the 
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magnetization direction Ena Hang. This approach is quite general and 
may be applied to express the relaxation law through the system parame¬ 
ters. In particular, it has been used to approximately describe the features 
of magnetic relaxation arising from the dipole-dipole interaction between 
particles Emma no] and rotating external magnetic field nn [ 22 ]. 

However, the described approach is rather complicated technically. The 
main difficulty is the necessity to take into account all possible directions 
of the magnetization. But the role of the states (magnetization directions) 
that differ appreciably from the equilibrium ones decreases with decreasing 
temperature. Therefore, if the total probability of these states is small 
enough, the state parameter of particles, which describes the relaxation 
process, can be approximated by the dichotomous random process. This 
so-called dichotomous approximation will permit us to study in detail the 
relaxation processes in a whole class of two-state systems. 

Within the dichotomous approximation, we associate the state param¬ 
eter of each structural element of the system with the dichotomous (or 
telegraph) process f(t), which takes the values —1 and 1 and satisfies the 
initial condition /(0) = 1, see Fig. [T] This process is characterized by wait¬ 
ing (residence) times Tj (j = l,oo), that is times between successive jumps 
of f(t), which are assumed to be independent random variables. We also 
assume that the waiting times t^-i in the up state, when f(t ) = 1, and 
the waiting times T^j in the down state, when f(t) = —1, are distributed 
according to arbitrary probability densities p + (r) and p _ (r), respectively. 



Fig. 1. (Color online) A sample path of the dichotomous process fit) with even 
number of jumps in the interval (0,f). 

The dichotomous process is the simplest random process that plays an 
important role in numerous applications. If the probability densities p ± (r) 
are exponential, then many of the properties of this process can be deter¬ 
mined by using the telegraph [if p + (r) = p~(r)] or generalized telegraph 
[if p + (r) / p~{t)\ equations [23, .Ml ES] ■ Since we are interested in study¬ 
ing the influence of different waiting time distributions on the character of 
relaxation of two-state systems, for this purpose it is more convenient to 
use the CTRW approach [23 27, Ml 122] ■ In this approach, we are going to 
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determine the probability density of the waiting times difference At, that is 
the difference between the times that the dichotomous process /(f) spends 
in the up and down states on the interval (0,t), defined as 

A t = dt' (2.1) 

Jo 

and study in detail the time dependence of the relaxation law 

MO = Pr{/(t) = 1} ~ Pr{/(t) = -1}, (2.2) 

where Pr{-} denotes the probability of the condition inside the braces. 

In order to avoid any confusion, we emphasize that all sample paths of 
the dichotomous process /(t) are assumed to start at f = 0 and /(0) = 1. If 
the sample paths started earlier, e.g., at f = —To with to > 0, then /(f) = 1 
as f € (—to,ti) and, due to ageing effects [30], the statistical properties of 
A t and the relaxation function /i(f) could strongly depend on to- Therefore, 
to simplify the problem, we restrict ourselves to the case To = 0. 

2.2. General results 

The probability density P(A, f) that the waiting times difference A* 
equals A for a fixed time f can be written as follows: 

P(A,f) = (<S(A t -A)>, (2.3) 

where <5(0 is the Dirac delta function and the angular brackets denote 
averaging over all sample paths of the dichotomous process /(f). Since 
Af € (—f,f), one gets P(A,f) = 0 as |A| > f, and the normalization condi¬ 
tion for P(A, f) becomes f_i t dAP(A , f) = 1. Let us assume that the process 
/(f) has exactly n (n = 0, oo) jumps in the interval (0, f). Then, introducing 
the n-jump probability density 

pW(A,f) = (i(A t -A))„ (2.4) 

((■) n denotes the average over these sample paths), we can represent P(A, f) 
in the form 

OO 

P(A,f) = ^pW(A,f). (2.5) 

72—0 

According to the definition (12.41) . the probability that a given sample 
path of /(f) has exactly n jumps in the time interval (0, f) is given by 

W {n \t) = f dApW(A,f), 


( 2 . 6 ) 
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and the normalization condition for the probability density P(A,t ) yields 
LoW^it) = 1. Our next step is to express W( n \t) and P( n \A,t) in 
terms of the waiting time probability densities p±(r). To this end, we first 
introduce the waiting times difference A^ for the sample paths with n 
jumps in the interval (0 ,t). It is obvious that Af 0) = t, and if n > 1 then 

n 

A<”> = + (-1)X +1 , (2.7) 

3 = 1 

where 

n 

T* +1 =t-J2 T j < T n+ 1- (2.8) 

3 = 1 

The probability that a sample path of the dichotomous process /(f) has 
no jumps on (0, t), that is the probability that t\ > t, is written as 

poo 

W^°)(f) = J drp + (r). (2-9) 

Let us now assume that the process /(f) has n > 1 jumps in the interval 
(0,f). If these jumps occur in the intervals (J2j=i T jil2jmi T j + drk) with 
k = l,n, then the probability dW^ n \t) of such process is given by 


dW( n \t) 




dTp n+ i(r), 


( 2 . 10 ) 


where Pj(r) = p + {r) or p~(r) if j is odd or even, respectively. It is appro¬ 
priate to recall here that Eq. ( 12 . 101 ) is obtained under conditions that (a) the 
jumps of /(f) are independent events with the probabilities pj(jj)dTj and 
(b) the (n + l)st jump occurs outside the interval ( 0 , f). To avoid confusion, 
we note that the times tj in Eq. ( 12 . 101 ) (and in all probabilistic expressions 
below) are interpreted not as the random variables but as the variables 
of integration. Introducing the re-dimensional domain of integration £l n (t) 
defined by the condition Ylj =i T j A t and replacing the lower limit of inte¬ 
gration in Eq. ( 12 . 101 ) by t* + 1 , we can write the probability W^ n \t) in the 
form (see also Refs. EH) 

wW(t) = / (YldTjPjiTj)) drp n+1 (T). ( 2 . 11 ) 

Finally, using the above results, we find 

P(A,f) = A) + P(A,f), 


( 2 . 12 ) 
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where P(A,t) = P^ n \A,t) is the regular part of the probability den¬ 

sity P(A,t) and 


P (n) {A,t)= [ (f[ dT jPj(' r j )][ d,Tp n+1 (r)6(A[ n) - A). (2.13) 

As it is seen from Eqs. (12.121) and (12.131) . P(A,t) depends on the waiting 
time densities p + (r) and p~(r) in a complicated manner. It seems reason¬ 
able therefore to find the connection between P(A,t) and p±(r) in Fourier- 
Laplace space. For this purpose, we first define the Fourier transform of a 
function ip( A) as T'l^A)} = <p k = J™ 00 dAe lkA ip(A) ( —oo < k < oo) and 
the Laplace transform of a function ip(t) as £{ijj(t)} = ij) s = / 0 °° dte~ st ip(t) 
(Res > 0). Then, taking the Fourier-Laplace transform of P(A,t) defined 
as Pks = C{F{P(A, t)}}, from Eq. (12.121) one obtains 

Pks = 1 ~ V _ s 7^ k + Pks , (2.14) 

where P ks = P^ and, according to Eq. (12.131) . 



To calculate P k ™\ we use the formula A = (— 1 ) n t — Y?j= i[(—l) n + 
(—1 yjrj, which follows from Eqs. (12.71) and (12.81) . and represent the inner 
integrals by the formula 


/•t—n. 


lQ n (t) j = \ 


~[dTjPj(Tj)= / dTipi(ri) / dT 2 p 2 (r 2 ) 




With these results, a straightforward integration in Eq. (12.151) yields 


dT n p n (T n ). 
(2.16) 


P, 


(2m—1) _ , + 


ks 


= (Ps-ik) m (Ps+ik) 

p(2m) _ / + - Ps-ik 

r ks ~ \P s—ikP s+ik) 


— 1 ^ P s+ik 

s + ik 


(2.17) 


s — ik 
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(m = l,oo). Finally, using Eqs. (12.171) and the formula for the sum of an 
infinite geometric series, rn = r /(l — r ) (M < 1), the Fourier-Laplace 

transform of P( A, t) can be written in the form 


Pks = 


Ps-ik 


( 


1 — p 


s-\-ik 


1 -Ps-ikPs+ik\ S + ik 


+ 


1 - pt-ik - \ 
s-ik Ps + ik J' 


(2.18) 


Note that Eq. (12.181) with p ± (r) = p(r) has been derived and used to de¬ 
termine the long-time behavior of P(A,t) in some particular cases (32j . 

Equations (12.141) and (12.181) represent the desired probability density 
P(A,t) in the Fourier-Laplace space. Because of the complex dependence 
of Pp s on k and s, the calculation of P(A,t) by taking the inverse Fourier- 
Laplace transform of Pk s is possible only in exceptional cases. In particular, 
if f(t) is the generalized telegraph process characterized by the exponential 
waiting time density functions 


p±( T ) = \±e~ x±T (2.19) 

(A± > 0 are rate parameters), then 

P( A, t) = e~ x+t 5(t - A) + ^ exp(- t - a) 

/A + A_(P — A 2 ) ) + A+A_(P - A 2 ))], (2.20) 


x 


where /o(') and P(') are the modified Bessel functions of the first kind of 
order zero and one, respectively. It should be noted that this result was 
previously obtained using other methods (see Refs. mmm)- 

We are also concerned with the time dependence of the relaxation func¬ 
tion p(t). According to the definition (12.21) . this function satisfies the con¬ 
ditions p(0) = 1 and p(t) € [—1,1]. Then, since (A t ) = J ( | dt'(f(t!)) and 
( f(t )) = p(t), it can be represented in the form p(t) = d(A t )/dt, which is 
useful if (A*) as a function of time is known. But to study the time depen¬ 
dence of fx(t) on waiting time distributions, it is more convenient to use the 
Laplace representation of /i(t). To derive it, we first express the relaxation 
function in terms of the probabilities W^ n \t)\ 


//(f) = IT (0) (f) + £ [W (2m) (t) - (t)]. (2.21) 

m =1 


Then, applying the Laplace transform to Eq. (12.211) and taking into account 
that wi°' ) = (1 — pf)/s and W^ n \t) = pj^ n \t )| fc=0 , we obtain 


-i + OO 

1 ~Ps + vA ,p(2m) _ p(2m-l) 
. ks ks 

m= 1 


k= O' 


Os = 


S 


( 2 . 22 ) 
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Finally, using Eqs. (I2.17p and performing summation over m , Eq. ([2.221) can 
be reduced to 


1 - 2 p+ +pfp s 
s(l -pfpj) 


(2.23) 


This result, which holds for arbitrary waiting time distributions, is our 
main tool for studying the relaxation of two-state systems. Before we pro¬ 
ceed to the analysis of relaxation laws, it is worthwhile to note that in the 
case of biased relaxation [when p + {r) / p“(r)] Eq. (12.23[) is equivalent to 
the integral equation 


p(t) — [ drp{r) [ cIt'p + (t')p (t — t — t') 

Jo Jo 

pt pt P~t — T 

= 1 — 2 d,Tp + (r)+ / (Jtp + {t) / (It'p~{t 1 ). (2.24) 

Jo Jo Jo 


In the case of unbiased relaxation [when p^{r) 
(12.241) are simplified to 


Ps = 


1 -Ps 
S(l+p s ) 


p(r)] Eqs. ([2.231) and 


(2.25) 


and 


p{t) + / dTp{r)p{t - t) = 1 - / dTp(r), 


(2.26) 


respectively. Equations ()2.24l) and (|2.26[) show that, in general, the memory 
effects in two-state systems play an important role for both biased and 
unbiased regimes of relaxation. 


3. Relaxation laws 

3.1. Exact results 

Here we consider the two-state systems characterized by Erlang distri¬ 
butions of waiting times, whose probability density functions are given by 


p ± (r) = 


X It 


k -1 


(fc-1)! 


= -A±r 


(3.1) 


where k = l,oo is the shape parameter. Since many of the properties of 
such systems have already been studied [35], we are only concerned with 
the relaxation function p(t). The conditions A_|_ 7 ^ and A_|_ = A_ = A 
correspond to the biased and unbiased relaxation, respectively. 
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3.1.1. Biased relaxation at k = 1 

At k = 1 the Erlang density functions (13.11) become exponential. Ac¬ 
cording to Eq. (12.191) . in this case pf = A±/(A± + s ) and Eq. (12.231) yields 


s — A + -)- A_ 

^ s s(s + A_)_ + A_) 


(3.2) 


By applying to Eq. (13.2|) the inverse Laplace transform [see Ref. [36], Eq. 5.2(5)], 
we make sure that the relaxation function is purely exponential: 


/i(f) 


A_ — A+ 

aTTaI 


2A + r ~(A + +A -)t 

\ + + X- 


(3.3) 


It can be easily verihed that this function satishes the differential equation 

—/r(t) + (A + + A ~)p(t) + A + - A_ = 0 (3.4) 

(/r(0) = 1), which also follows from Eq. ([2.24D . Thus, there are no memory 
effects in this case. 


3.1.2. Biased relaxation at k = 2 

At k = 2 from Eqs. (13.11) and (12.231) one can get pf = A±/(s + A±) 2 and 


s 3 + 4as 2 + 4(i^ 2 + ab)s + 4 u 2 b 
s(s + S-|_)(s + s_)(s + 2a) 


where a = (A+ + A_)/2, b = (A_ — A+)/2, v = A_|_A_, and s± = a ± 

\Ja 2 — 2 v 2 . Taking the inverse Laplace transform of Eq. (13.51) [see Ref. [36 ]. 
Eq. 5.2(19)], for a given class of two-state systems we obtain 


p(t) 


X - ~ A + (i ^± P -(A++A_)t 

A + + A_ 1 + A_ 6 



cosh (va 2 


2v 2 1 ) 


— (A + — 3A_ 


sinh (Va 2 — 2v 2 t)' 
2 s/a 2 - 2v 2 . 


g—(A++A- )t /2 


(3.6) 


Although the relaxation functions (13.31) and (13.61) have the same limiting 
value n{ oo) = (A_ — A + )/(A + + A_), their behavior at finite times is quite 
different. In particular, at small times the function l — p(t) is proportional to 
t and t 2 , respectively. Moreover, in contrast to (13.31) . the relaxation function 
(13.61) at a 2 — 2v 2 > 0 is characterized by three relaxation times, the largest 
of which is l/s_. 
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The most important qualitative difference between the relaxation laws 
m and (13.61) occurs if a 2 — 2v 2 < 0, that is if the parameters \± satisfy 
the conditions 

3 - 2y/2 < < 3 + 2V2. (3.7) 

A_ 

Since in this case y/ a 2 — 2u 2 = i\J2v 2 — a 2 with 2v 2 — a 2 > 0, the hyper¬ 
bolic functions in Eq. (13.61) should be replaced by the trigonometric ones. 
Therefore, introducing the period of these functions as T = 2-k/ y/2v 2 — a 2 
or 


T = 


47T 

^/6A+A_ - A 2 . - A 2 ’ 


(3.8) 


Eq. (13.61) can be rewritten in the form 




cos 


/ 27rt 

\T~ 


,, , , T (2nt \ 

-(A + -3A_)— sin (—) 


s —(A++A_)i/2 


(3.9) 


According to this formula, the function p(t) tends to the limiting value 
p(oo) in an oscillating manner (if the conditions (13.71) hold). This is a 
remarkable and somewhat unexpected result because the waiting time den¬ 
sities do not contain periodic functions. In view of our previous discussion, 
such behavior of /r(t) can be interpreted as the emergence of memory ef¬ 
fects. The non-monotonic dependence of /r(t) on t obtained from Eq. (13.91) 
is illustrated in Fig. [2] for different values of the rate parameters X±. 


3.1.3. Unbiased relaxation at k > 2 

In this case p s = A k /(s + \) k and, according to Eq. (12.25[i . we obtain 


1 — (s + X) k /X k 

[1 - (s + A)/A] [1 + (s + X) k /X k ] • 1 j 

The last result shows that p(t) = e~ xt g(Xt ), where the function g(t ) is 
defined by its Laplace transform 


9 s = 


1 - s k 

(1 — s)(l + s k ) ’ 


(3.11) 


Using the relations (1 — s k )/( 1 — s) = Ylj=i 1 an d 1 + s k = n^=i(« — «?) 
with aj = e l ' n ( 2 i~ l )/ k being the solution of the equation 1 + s k = 0, we can 
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t 


Fig. 2. (Color online) The relaxation functions in the case of Erlang distributions of 
waiting times for k = 2, A + = 0.5, A_ = 1 (1), A + = A_ = 1 (2), and A+ = 1, A_ = 
0.5 (3). The solid lines represent the theoretical result (13.911 and the symbols 
show the simulation results for ( f(t )). The period T of trigonometric functions in 
Eq. (13.911 equals 87 t/-\/7 and 27r for the cases (1), (3) and (2), respectively. 


rewrite the above formula as 


EjU s j ~ 1 

n ; ; ,(* 


(3.12) 


Then, taking the inverse Laplace transform of Eq. (13.121) [see again Ref. [35] , 
Eq. 5.2(19)], we find 


9(t) = E 


£*=i of 


j-i 


,<Ht 


1=1 n j=i ( a i a j ) 


(3.13) 


where the prime on the product means that j ^ l. This result, together 
with the fact that J2j =l c^ -1 = 2/(1 — a;), yields 






i=i 


e -(l -ai)Xt 

(l-o,)n?=i(ai-ai)‘ 


(3.14) 


Finally, since nj=i( a z — a j) = kaf 1 = ~k/ai, the relaxation function (13.1411 
can be represented in the form 


/r(f) = 


Ok - 


2\t 


[k/2] r 


+ iE 


1=1 


cos (A At) 


A 


1 - Oil 


sin(AAf) 




(3.15) 


Here, Ok = 
of k/2 , and 


0 or 1 if is even or odd, respectively, [k/2] is the integer part 


Oil = cos 


(—*—)’ = <3 - 16) 
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Thus, according to Eq. (13.151) . the relaxation to the final state [p(oo) = 0] 
also occurs in an oscillating way. But, in contrast to the relaxation function 
(13.91) that oscillates with a single period T, in this case the oscillating part of 
p(t) is, in general, characterized by a few periods T) = 27r/(/3;A) (if k > 5). 


3.2. Asymptotic results 

Next we study the long-time behavior of the relaxation function in two 
types of systems. One of them is characterized by heavy tails of waiting 
time distributions, and the other by superheavy tails. Our aim here is to 
find universal asymptotic laws for unbiased relaxation, when p±{r) = p(r). 


3.2.1. Heavy-tailed p(r) 

A class of heavy-tailed probability density functions of waiting times is 
defined by the asymptotic behavior 

P( r ) ~ 7ITT ( 3 ‘ 17 ) 

(t —>• 00 ), where q is a positive parameter and the tail index a satisfies the 
condition a € (0,2]. Since the long-time behavior of p(t) is related to the 
behavior of p s when the real parameter s tends to zero [371 [38], we need to 
find p s as s —>• 0. Using the asymptotic formula (13.171) . it is not difficult to 
show (see, e.g., Ref. [39]) that 




1 - Ps ~ < 


qs In i 
ts — q 

. ts - fs 2 hii, 


r( 2 — a) na 

TS ~ 9 Q (a-l) S , 


OL G (0, 1) 

a = 1 
a € (1,2) 
a = 2 


(3.18) 


(s —y 0), where T(-) is the gamma function and f = / 0 °° drrp(r). Based on 
this result, the asymptotic behavior of pit) can be obtained from the Taube- 
rian theorem for Laplace transforms. It states [371 [38] that if a function h(t) 
is ultimately monotone and 


h s ~ 



(3.19) 




as s —>• 0 then 


h{t) 


r^j 


(3.20) 
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as t —>■ oo. Here, p > 0 and L{t) is a positive function slowly varying at 
infinity, i.e., function for which the condition L(at) ~ L(t) (t —>• oo) holds 
for all a > 0. 

If ct € (0,1) then, according to Eqs. (12.2511 and ([3.181) . p 3 ~ gT(l — 
a)/(2as 1_ “) (s —> 0) and, associating /i(f) with p(t), directly from the 
Tauberian theorem one gets 


(3 ' 21) 

as t —>• oo. Since at a = 1 and s —>• 0 the condition // s (g/ 2 )ln(l/s) 
occurs, the Tauberian theorem is not applicable to p(t). However, it can 
be applied to the function h(t) = Jj dt'p(t') whose Laplace transform is 
given by h s = p s /s. This yields h(t) ~ (q/2)lnt as t —>• oo and, because 
//(f) = dh(t)/dt, one can make sure that the case a = 1 is described by 
the asymptotic formula (13.211) as well. Similarly, introducing the auxiliary 
function h(t) = r /2 — Jq dt'p^t') and using the relation //(f) = —dh(t)/dt, 
it is not difficult to verify that the above formula holds also for a € ( 1 , 2 ]. 
Thus, any two-state system with a heavy-tailed distribution of waiting times 
exhibits slow (power-law) relaxation (13.211) . 

3.2.2. Superheavy-tailed p(r) 

Let us finally consider a class of superheavy-tailed probability density 
functions of waiting times defined by the asymptotic behavior 

(3.22) 

T 

(t —>• oo). Here, Z(r) is a slowly varying function which, due to the nor¬ 
malization condition f£° drp{r) = 1 , tends to zero in such a way that 
1(t) = o(l/lnr) as r ^ oo. This class of probability densities has re¬ 
cently been used to study the phenomenon of superslow diffusion within 
the CTRW formalism |4TI. ;il, 121133] , 

To find the long-time behavior of the relaxation function p(t) in such 
two-state systems, we first introduce the exceedance probability V(t) = 
f t °° drp(r). Its important feature, which follows directly from Eq. (13.221) . 
is that it tends to zero at t —» oo as a slowly varying function. This fact 
enables us to determine a universal asymptotic formula for p{t). Indeed, 
since 1 — p s = sV s and V s ~ (l/s)V(l/s) (s —> 0), Eq. (I2.25|) yields p s ~ 
(l/2s)V(l/s). From this, using the Tauberian theorem, one obtains 

Kt) ~ \ V( J) 


(3.23) 
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as t —>• oo. Since the condition lim^oo t p V(t) = oo holds for all p > 0 [33] . 
p(t) decreases slower than any negative power of time. It is therefore the 
relaxation process, whose long-time behavior is described by Eq. (13.231) . can 
be called superslow. In particular, if 


p{r) 


Inc 

(c + t) In 2 (c + t ) 


(3.24) 


(c > 1) then the exceedance probability equals V(t) = In c/ ln(c + t) and 
Eq. (13.231) gives 

. . Inc , 

m ~ (3.25) 

The asymptotic behavior of the relaxation function p(t) in the cases of 
heavy and superheavy tails of the waiting time distribution p(r) is illustrated 
in Fig. [2] It shows that, because V(t) is a slowly varying function, the 
function p,(t) in the regime of superslow relaxation tends to zero much slower 
than that in the regime of slow relaxation. 



Fig. 3. (Color online) Asymptotic behavior of the relaxation function p(t) in 
regimes of slow (1) and superslow (2) relaxation. The line 1 represents the asymp¬ 
totic result (13.211) with <7=1 and a = 1, and the line 2 represents the asymptotic 
result (13.251) with c = 2. 


4. Conclusions 

Using the continuous-time random walk approach, we have studied the 
phenomenon of relaxation in a class of two-state systems, whose structural 
elements evolve according to the dichotomous process. Our interest has 
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been focused on the probability density function of the waiting times differ¬ 
ence and on the relaxation law. Assuming that the distributions of waiting 
times in the up and down states of the dichotomous process are arbitrary, 
we have found the Fourier-Laplace representation of this density function 
and the Laplace representation of the relaxation function. These represen¬ 
tations have been used to determine the density function in the case of 
exponential distributions of waiting times, to derive the integral equations 
describing the biased and unbiased relaxation processes, and to calculate the 
relaxation laws in some special cases. In particular, we have considered the 
Erlang distributions of waiting times and two classes of distributions charac¬ 
terized by heavy and superheavy tails, respectively. It has been shown that, 
depending on the waiting time distributions, the two-state systems can ex¬ 
hibit a wide variety of relaxation regimes including exponential, oscillatory, 
slow and super slow ones. 
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